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The IMS -H I S double-focussing mass spectrometer that flew on the 
GIOTTO spacecraft covered the mass per charge range from 12 to 56 
(AMU/e). By comparing flight data, calibration data and results of 
model calculations of the ion population in the inner coma, the 
absolute mass scale is established, and ions in the mass range 
25 to 35 are identified. Ions resulting from protonation of mole- 
cules with high proton affinity are relatively abundant, enabling 
us to estimate relative source strengths for H 2 C0, CH 3 OH, HCN, and 
H 2 S, providing for the first time a positive in situ measurement of 
methanol. Also upper limits for NO and some hydrocarbons are derived. 
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The composition of gases in the coma of comets can be studied by 
spectroscopy at various wave lengths of electromagnetic radiation, by 
mass spectrometric analysis of the neutral component, and by ion mass 
spectrometry. All these methods have their advantages and limitations 
and it is their combination that has led to the advances in cometary 
chemistry that we have witnessed since the flight of spacecraft into 
the coma of Halley's comet. 


In this paper, we use the data obtained by the HIS sensor of the 
IMS experiment (Balsiger et al„ 1986) on the spacecraft GIOTTO during 
its passage through the inner coma of Halley's comet. The HIS is a 


double-focussing magnetic mass spectrometer with 9 channel electron 
multipliers (OEMs) arranged in a specially designed distributor, 
called an “Igel", at its exit (cf. Balsiger et al., 1986). It provides 
information on the mass/charge ratio, the velocity and the angles of 
incidence of the ions. Among the GIOTTO instruments, the capabilities 
of the HIS are unique for investigating the ion composition 
region between the cometopause (Gringauz et al., 1986) and the lono- 
pause (Neubauer et al., 1986; Neubauer, 1987; Balsiger et al.. 1986a) 

of Halley's coma. 


The HIS, of course, also works in the low temperature (200-300 K) 
region inside the ionopause, but it is not simple to reduce the four- 
dimensional representation of these data (i.e. mass/charge, velocity 
and two angles of incidence) to a one-dimensional mass spectrum giving 
quantitative ion abundances, because the spectrometer response is not 
strictly a separable function of these four parameters. Thus, the 
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apparent mass/charge ratio depends somewhat on the velocity and direc- 
tion of incidence. We have developed a method by which quantitative 
M/Q spectra are derived, using a scheme of adding CEM count rates 
combined with a fitting routine based on laboratory calibration data 
(Meier, 1988). 

In contrast to previous attempts to reproduce the IMS mass 
spectra by ion-chemical modelling (e.g. Wegmann et al., 1987; cf. 
Schmidt et al., 1988), we use a model which concentrates on the ion 
mass range between 25 and 35 AMU/e in order to obtain better esti- 
mates of the parent molecule abundances which contribute to this ion 
mass range. 


2. Experimental Results 

The HIS sensor took data in two modes, the H- and the N-mode (Balsiger 
et al., 1986; Meier, 1988). Results obtained in the H-mode for the 
mass range 12 to 35 at nine distances from the nucleus are shown in 
Figs. 1 and 2. Each spectrum represents the data of one spin period, 
corresponding to a traversal of 273 km. The distances from the nucleus 
given in the figures are determined at the centers of the respective 
spin periods. The HIS instrument scans the voltage of the energy 
analyzer In 64 steps, covering the energy range appropriate for rammed 
cold ions with mass/charge 12 to 56. In the innermost part of the 
coma a countrate vs. measurement step plot represents a mass spectrum, 
because all cold ions have the same velocity in the instrument frame 
of reference, i.e. the ram velocity (68.4 km/s). The H- and N-mode 
are different in so far that in the H-mode all masses with the same 
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angle and velocity are measured by the same CEM, whereas in the N-mode 
different masses are deflected to different CEMs. Hence the H-mode is 
well suited to get full survey mass spectra (by plotting countrates 
vs. measurement step). However, the assignment of mass in the measure- 
ment step/CEM-matrix is dependent on ion speed, temperature and angle, 
and the detailed response of the matrix bins to each ion beam had to 
be carefully calibrated. Thus, for the purposes of this paper, we use 
the H-mode to give an overview of the spectra as a function of dis- 
tance from the nucleus. The N-mode data are then used to derive rela- 
tive ion abundances based on the calibrated response function. 

Up to about mass 20, the identification of the ionic masses does 
not pose a problem, thus there is no uncertainty in the mass scale 
given in Fig. 1. This is supported by the fact that the calibration 
data predict the position of HsO + where the highest ion peak is actual- 
ly found. From this part of the spectrum, we can determine velocity 
and angle of incidence of the ions. If we then assume that the velo- 
cities and the flow directions are the same for all masses, we can 
assign an accurate mass scale to the higher masses as well. In Fig. 

2, we give the H-mode spectra thus obtained for the ions in the range 
M/Q = 25-35. 

The data displayed in Figs. 1 and 2 divide naturally into three 
spatial regions as can be recognized in the patterns of the spectra. 

I. 1215 km - 3289 km . This region lies well inside the contact 
surface (Neubauer et al., 1986; Neubauer, 1987; Balsiger et al., 

1986a). Here the counts are fairly evenly distributed in azimuth 
indicating that the flow direction coincides with the spin axis. The 
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ion temperature is less than 500 K. The M/Q patterns of the five 
spectra are fairly similar and indicate a good tuning of the in- 
strument and a good mass resolution for the regular flow velocity 
encountered here. 

II. 4365 - 6554 km . This region lies mainly outside the contact 
surface (at 4660 km, Neubauer, 1987). The counts are unevenly distri- 
buted over azimuth in this region, indicating a flow with a velocity 
component coming from the solar direction. This is consistent with 
expectation since bulk flow of the ions should follow roughly the 
contours of the contact surface. The 4365 km spectrum contains data 
from inside the contact surface. It was here that the HERS sensor 
detected hot cometary ions wjiich may be due to the double charge- 
exchange mechanism (Goldstein et al., 1987; Eviatar et al., 1989). 

The flow conditions and distribution functions in region II are sig- 
nificantly different from those in region I. The poor mass definition 
and resolution shown for the region II spectra in Figs. 1 and 2 
results from the fact that they were generated using the flow con- 
ditions prevailing in region I. 

III. 9814 km . This spectrum is typical of the region outside 
^7500 km. The azimuthal anisotropy is still observed, compatible with 
the expected flow pattern around the contact surface. The spectrum is 
similar to those in region I, with the exception of the strong count 
rate at M/Q * 32. The latter may result from a change in chemical 
composition with distance (i.e. there is a slow build-up of S + ; cf. 
section 4). 



3. Interpretation of the spectra 
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This paper is concerned primarily with the ion mass range 25 to 35 
AMU/e (Fig. 2) focussing on region I, i.e. the low temperature gas 
inside the contact surface. Here, photodissociation, photoionization, 
ion-molecule reactions, and dissociative recombination dominate the 
chemistry. The dominance of water vapor in the neutral gas together 
with its high proton affinity (PA) results in H30 + being by far the 
most abundant ion. Thus, this ion governs the direction that the ion- 
molecule reactions take (Aikin, 1974). As a consequence, the PA of 
the molecules or radicals in the inner coma (Table 1) is the most 
Important chemical parameter (Huntress, 1977), determining what the 
main ionized species are for a particular parent molecule. 


Molecules (M) having a higher PA than water (cf. Table 1) react 
readily with H 30 + abstracting a proton: 

H 3 0 + + M -*■ H 2 0 + MH + . (kl) 

Consequently, the resulting MH + ions (Table 2, column 3) are very 
abundant relative to their parent molecules. For molecules containing 
H, C, and 0, reactions with HsO + have the effect of transforming a 
molecular mass spectrum dominated by even numbers into an ion mass 
spectrum dominated by odd numbers. 

On the other hand, ions M + produced by simple photoionization 
(emission of an electron) from a molecule having a PA lower than 
water fall into two categories. If they react with water (e.g. C0 + , 
cf. Table 1) 


M + + H 2 0 - 


(k2) 


products 
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they will have a low abundance (Table 2, column 5). If they do not 
react with water (e.g. NO* ; cf. Table 1), they will have a rather 
long lifetime, governed in most cases by dissociative 
and thus they will accumulate towards a relatively high abundance 

(Table 2, column 4). 

While these simple systematics are qualitatively valid, 
molecule reactions with species other than the most abundant H 2 0 and 
H 3 o + must be taken into account for a quantitative analysis. For this 
purpose, we have developed a numerical method which is simple, yet 
adequate for treating the radial evolution of ion abundances in the 
inner coma (Appendix A). The 21 neutral and 40 ion species listed in 
Table 3 were included in the calculations for the present paper. The 
parent molecules selected fall into three categories. (1) Those whic 
could reasonably be expected and are essential to obtain an adequat 
fit to the data (e.g. H 2 0, CO), (2) those which have been previously 
reported but cannot be positively confirmed or excluded by the data 
(e.g. N 2 ), and (3) theoretically possible molecules for which the 
dau provide abundances (H 2 C0, CH 3 0H, NON, H 2 S) or meaningful upper 
limits (e.g. NO, hydrocarbons). Photodissociation and photoionization 
of all neutral species, dissociative recombination of the ions, and 
about 280 ion-molecule reactions were taken into account. Information 
on the rate constants that were used is given In Appendix B. While 
the expansion velocity could be introduced into the program as a func- 
tion of radial distance, for the purpose of this paper, which deals 
with the coma Inside the contact surface, the expansion velocity was 
assumed to be a constant 900 m/s (cf. Lammerzahl et al., 1587). The 
program allows for molecular point sources (near the nucleus) and for 
extended sources. We have used a point source for B 2 0 with a strength 
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corresponding to a H 2 0 density of 5*10 7 cm 3 at 1000 km from the 
nucleus (Krankowsky et al., 1986). Strengths of point sources or ex- 
tended sources for the other neutral species given in Table 3 could 
be freely chosen. Many numerical integrations varying these source 
strengths were carried out in order to study the influence of indi- 
vidual components and combinations of components. 

As discussed above, the mass scale can be influenced by various 
parameters including angle of incidence and ion speed. The results 
can be summarized as follows: With the mass scale given in Fig. 2, 

we find a reasonable set of molecular abundances that can reproduce 
the observed ion mass spectra in the mass per charge range 25 to 35. 
We have not been able to find an alternative mass assignment that 
could do this. For instance, a shift of all peaks by one unit to 
higher M/Q values could not be ruled out absolutely by the calibra- 
tion data alone. However, we would have no explanation for a high 
abundance at M/Q = 34 along with a low abundance at M/Q = 35; e.g. 
H 2 S would produce the opposite effect. Also, a relatively high flux 
of M/Q = 29 ions could not be explained: HC0 + has to be lower than 

C0 + , and C 2 H 5 + or N 2 H + react with H 2 0 with high rate constants. 

A shift of all peaks by one unit to lower M/Q values would be 
difficult to reconcile with the calibration data. Moreover, it would 
again be difficult to find a chemical interpretation (cf. Table 2). 
For instance, it would be hard to explain the high abundance at 
M/Q = 30 that results with this mass scale. 

On the other hand, the mass scale given in Fig. 2 results in 
relative abundance peaks inside the contact surface at the four mass 
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numbers corresponding to the products of H 3 0 + listed in column 3 of 
Table 2. In other words, the pattern in the mass range 25 to 35 finds 
a natural explanation. A careful comparison of this mass scale with 
the calibration data implies that the average arrival direction of 
the ions was M 3. 3 degrees away from the optical axis of the entrance 
deflector of the HIS sensor and that the ions had an average radial 
speed of 900 m/s. These parameters also fit the lower mass ions given 
in Fig. 1, for which the mass assignment is not in question. 

With the adopted mass scale, we interpret the spectra obtained 
between 1215 km and 3289 km as follows: 

1. Because of their high proton affinity, H 2 CO and CH 3 OH dominate 
the mass range 25 to 35, being the parent molecules of the two most 
abundant ions (CH 3 0 + and CH 3 0H 2 + ). 

2. The minimal ion abundances at masses 27, 29, and 34 are con- 
sistent with prediction. All ions listed in Tables 2 and 3 with these 
masses (C 2 H 3 + , HCN + , C 2 H 5 + , N 2 H + , HC0 + , H 2 S + ), and also H 2 0 2 + , react 
with water. 

3. The peak height at mass 28 is consistent with the HCN/H 2 0 
abundance ratio of 0.001 derived by Schloerb et al. (1987) from radio 
frequency observations. C0 + also contributes to the mass 28 peak but 
is not dominant. 

4. The hydrocarbons C 2 H n are difficult to determine from the ion 
spectra. At this time, only upper limits can be given. The small peak 
at mass 26, if real, would be due to C 2 H 2 + , rather than CN , since 
the latter is readily destroyed by water (Tables 1 and 2). 
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For the purpose of deri ving-molecular abundances, we have com- 
puted ion mass spectra as a function of distance from the nucleus, 
varying the proportions of the assumed parent molecules to optimize 
the fit. In Fig. 3, we compare such a computed spectrum, based on the 
abundances indicated, with the observed spectra at two radial dis- 
tances straddling that of the computation. It is evident that the fit 
is quite reasonable. The sensitivity of the fit to variations in the 
parent molecule source strengths was investigated in order to estab- 
1 ish ranges or 1 imits . 

As noted in Fig. 3, contributions from extended sources were 
assumed for some molecules. For these we assume an exponential de- 
pletion length of 10 4 km, originating at the nucleus. Such extended 
sources are indicated by the comprehensive CN-observations of A'Hearn 
et al. (1986) and Hoban et al. (1988), as well as by the mass spectro- 
meter results for neutral CO (Eberhardt et al., 1987). While the need 
to include both point and extended sources is clearly indicated, both 
from the previous work and our data, the distance between the two 
experimental spectra given in Fig. 3 is too small for arriving at a 
reliable quantitative separation between them. For this, a complete 
Investigation of the radial evolution of each ion in question would 
be required. 

We emphasize that the molecular abundances shown in Fig. 3 are 
not unique in giving a reasonable agreement between theoretical and 
experimental spectra, considering the uncertainties in the experi- 
mental data and in the reaction rate constants. Definite abundance 
estimates can be obtained, however, for molecules having high PAs, 
because the ions resulting from their protonation dominate the mass 
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spectra. This is evident from Table 4 where the relative contribu- 
tions to some relevant mass numbers are listed. 


4. Molecular Abundances 


Below, we give abundance estimates derived from a comparison of the 
results of the model calculations with the measured spectra in region 
I. We do not distinguish between point sources and extended sources 
for most molecular species. Some of the abundances might be improved 
by applying a more refined fitting programme to the data and by com- 
paring experimental and theoretical ion abundances over the whole 
length of the trajectory inside the contact surface. However, we are 
confident that the resultant improvements would not materially alter 
the conclusions, and thus the added complexity of the modeling is not 
justified for the purposes of this paper. 

Ammonia. NH 3 can be estimated from the NH 4 /H 3 O ratio (Ip, 
1986). Allen et al. (1987), using IMS data, applied this method to 
derive a NH 3 /H 2 O abundance ratio of .01 to .02; however, they neg- 
lected reactions of H 30 + with molecules heavier than water. Ip et al . 
(1990) pointed out that even relatively rare HCN influences the 
ammonia abundance thus derived. Now that we find a few to several 
percent of molecules with higher PA than water, we obtain a reduction 
of H 3 0 + , which in turn somewhat reduces the estimate on the NH3/H2O 
ratio. The question of whether or not the ratio thus obtained is com- 
patible with the NHs/HaO value of 0.003 derived from the NH 2 emission 
(Wyckoff et al., 1988) will be discussed in a separate paper. 
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Formaldehyde. The mass 31 ion rate falls off less rapidly than 
R ” 1 between 1215 km and 9814 km (cf. Fig. 2) in spite of the relatively 
short photo lifetime of H 2 CO (cf. Table 2). From this, we derive a H 2 CO 
source which is at least partly extended with an integrated strength 
(relative to water) of about one percent at 1500 km, rising to a few 
percent at the ionopause. The mass spectra measured in region II sug- 
gest that the source extends even farther out. Our results are in 
agreement with data obtained by other instruments. H 2 CO has been found 
with the IKS infrared spectrometer on the VEGA 1 spacecraft (Combes 
et al., 1988; Mumma and Reuter, 1989). Polymeric formaldehyde has 
been inferred from the spectra obtained by the ion energy spectro- 
meter PICCA onboard the GIOTTO spacecraft (Mitchell et al, 1987; 
Huebner, 1987) . Krankowsky (1990) found the formaldehyde molecule 
with the GIOTTO NMS mass spectrometer. He estimated a H 2 CO abundance 
relative to water of 4.5 percent, including a contribution from an 
extended source. 

Methanol. The peak at mass 33 could be protonated methanol or 
hydrazine. Stief and deCarlo (1965) and Delsemme (1975) have proposed 
that the 3360 AngstrSm emission observed in comets is due to a NH 
radical, the precursor of which is hydrazine. Linder the thermodynamic 
conditions In the solar nebula or the proposed sub-nebula (Fegley and 
Prinn, 1988), N 2 H 4 would not be produced in significant amounts. Cal- 
culations show that the equilibrium ratio N 2 H 4 /N 2 for these two sites 
is only 2*10‘ 38 and 3*10' 29 at 300 K or 4*10' 20 and 3*10" n at 1500 K. 
The equilibrium ratio CH 3 OH/CO is higher, especially in a sub-nebula 
environment. Also, there has been no mentioning of hydrazine produc- 
tion in interstellar clouds, whereas some ammonia and methanol are 
produced according to ion molecule reaction calculations (cf. Prasad 
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and Huntress, 1980; Leung et al., 1984). Moreover, if we wanted to 
identify the mass 33 peak with N 2 H 5 + , we would have to postulate 
N 2 H 14 /NH 3 £ 1 which seems implausible. Thus, we propose that the mass 
33 peak is essentially due to CH 30 H 2 + which leads to an estimate of 
the methanol/ water ratio of .003 to .015. This is the first in situ 
measurement of CH 3 0H in a comet. However, methanol has also recently 
been found in comet Austin by millimeter observations (Colom et al., 
1990). 


Nitrous Oxide. The count rate at mass 30 places an upper limit 
of .005 on N0/H 2 0. NO is usually not mentioned among possible cometary 
constituents. Since, however, theory indicates that NO is produced in 
dense interstellar clouds (Prasad and Huntress, 1980), it may be well 
to mention this NO abundance limit for Halley's coma. 

Hydrogen Cyanide. The peak at mass 28 is probably mainly due to 
H 2 CN + . Although CO is the second most abundant molecule in Halley's 
coma (Eberhardt et al., 1987), the contribution of C0 + is minor, be- 
cause this ion reacts readily with water (cf. Tables 1, 2, and 4). If 
we assume that also C 2 H 4 does not contribute very much, then we have 
enough H 2 CN + to account for the HCN/H 2 0 ^ .001 abundance ratio derived 
from radio frequency observations (Schloerb et al., 1987) or from the 
CN emission (A'Hearn et al., 1986). Our estimate for the HCN abundance 
is higher than the one given by Ip et al. (1990). 

Hydrocarbons. Since C 2 Hg, C 2 Hi,, and C 2 H 2 have PAs lower than 
water, they do not form ions with prominent abundances. However, the 
count rates at masses 25 and 26 indicate that some hydrocarbon ions 
are present in the mass range considered. While a careful analysis of 
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the abundances - particularly at even mass number - might produce some 
definite hydrocarbon estimates, we simply mention that as a result of 
cur present modelling effort the sum of the abundances of acetylene 
and ethylene are at most one or two percent relative to water. 

h 2 s and s. The inclusion of S, H 2 S and related ions was intended 
to allow for studying their contribution in the mass range considered. 
We defer a more complete discussion of the sulfur containing species 
to a later publication and discuss here only the abundance of H 2 S, 
which has just been detected in comet Austin by millimeter observation 
(Colom et al., 1990). For comet Halley's coma, Marconi et al. (1990) 
derived a lower limit of 0.5 percent for the H 2 S/H 2 O ratio. However, 
their conclusion is based on a dissociation time for H 2 S of ^300 sec, 
which we consider to be erroneous for two reasons: (1) The photolysis 
of H 2 S leads to SH which in turn dissociates in 118 sec (0.9 AU; Kim 
and A'Hearn, 1990) to give S. From atomic sulfur, S + is produced at a 
relatively high rate (Table Bl). Thus, not only H 3 S but also S + carries 
significant information on the occurrence of H 2 S. Our model calcula- 
tions clearly show that H 2 S with a half percent abundance and a life- 
time of only ^300 sec would lead to a much higher ion abundance at 
mass 32 than we actually observe (cf. Fig. 3). (2) Using the cross 
section data given by Lee et al. (1967), Kim and A'Hearn (1990) have 
derived a photodissociation time for H 2 S of 3200 sec (0.9 AU), a result 
we have confirmed. With this longer dissociation time, we have consis- 
tency between the masses 35 and 32 in the ion spectra and obtain a 
H 2 S source strength (point source and/or extended source) of 0.001 to 
0.004 relative to H 2 O. 
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There is a significant change in the appearance of ion mass 
spectra in region II relative to region I (Figs. 1 and 2). As pointed 
out in Section 2, the distribution functions of the ions have changed 
significantly from the quiet flow encountered in region I. This makes 
it difficult to identify changes in ion abundances with the methods 
used in this paper, since these methods are based on the assumption 
of quiet flow conditions appropriate to region I. 

While the perturbations in the spectra observed in the mass range 
16-20 in region II have again disappeared in region III (Fig. 1), the 
high count rate at M/Q = 32 remains (Fig. 2). A slow build-up of 
M/Q = 32 actually begins between 2000 and 3000 km. Thus, this general 
increase in M/Q = 32 with increasing distance could be due to the $ + 
ion. However, quantitative interpretation of the spectra in regions 
II and III needs further analyses, which we leave to future publica- 
tions. 


5. Discussion 


The ions CH 3 0 + and CH 3 0H 2 + are derived from formaldehyde and methanol, 
respectively. Judging from ion-molecule reaction systematics, it seems 
improbable that heavier gaseous aldehydes, ketones or alcohols contri- 
bute significantly to these ions. Thus, we propose that both H 2 C0 and 
CH 3 0H are present in the gas phase with the relative abundances indi- 
cated in Section 4. This observation neither precludes nor supports 
the possibility that H 2 C0 and/or CH 3 0H are present in some poli- 
merized form in the solid (Vanysek and Wikramasinghe, 1975; Huebner, 
1987; Mitchell et al., 1987). The radial decrease of CH 3 0 + is less 



17 


than expected for a point source, indicating that H 2 C0 is released 
from grains over times of the order of hours. 

The analysis offered here underlines that carbon in Halley's 
comet exists largely in oxidized form - at least the portion of it 
that is released into the gas phase within the first one or two hours. 

Krankowsky et al. (1986a) have shown that the dominating peaks 
in the 40-50 mass region occur at M/Q = 45 and M/Q = 47, and they 
identify them as HCS + and H 3 CS + , respectively. We note that protonized 
acetaldehyde (CH 3 CHO) and ethylalcohol (C 2 H 5 OH) may Significantly 
contribute to these mass peaks. 
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Appendix A 

The chemistry in the ionosphere - the region well inside the contact 
surface - of Halley's comet is dominated by photodissociation, photo- 
ionization, ion-molecule reactions, and dissociative recombination. 
Since these reactions and in particular the ion-molecule reactions 
connect a great variety of species in a rather non-systematic way, a 
large matrix of non-linear differential equations results. 

The ionosphere is characterized by a steady expansion velocity 
of ^-0.9 km/s which is roughly equal for ions and neutrals (Krankowsky 
et al., 1986; Lammerzahl et al., 1987) and by low kinetic temperatures 
(200-300 K; Balsiger et al., 1986a). Thus, the rate constants for 
ion- molecule reactions measured in the laboratory can be used. 

We found reactions between neutrals (molecules or radicals) to 
be of minor importance for the species considered here; thus, we have 
neglected this type of reaction. In fact, the most important change 
in neutral abundances comes from photodissociation (and Ionization) 
and this has been taken into account. 

The ion abundances originally obtained from photoionization are 
thoroughly changed by ion-molecule reactions. On the other hand, for 
the neutral partners the relative importance of these ractions is 
small because of the relatively low abundance of the ions. Thus, we 
have neglected the changes in neutral abundances due to ion-molecule 
reactions. The types of reactions included in our model are given in 
the following equations (reactants or products for which the changes 
due to the reaction are taken into account in our calculations are 
underl inded) : 
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A +hv + B + C + (1) 

A+hv-*-B + + e2+C+ (2) 

A + + C + (3) 

-*■ C + D (4) 


It is seen that with the simplifications introduced here, the 
first three types of equations are linear. 

The set of differential equations corresponding to the reaction 
equations (1) to (4) is solved by stepwise integration over r, the 
distance from the nucleus. The relation between time t and r Is given 
by the expansion velocity v(r) which was taken to be constant at 
900 m/s. A backward integration scheme was chosen for the destruction 
terms, but not for the production terms which were treated according 
to the forward method. In this way, the differential equations are 
decoupled, and a matrix conversion procedure is avoided. 

The reaction rates of equations (1) to (4) form one-, two- or 
three-dimensional arrays of real numbers. Since many elements in 
these arrays are empty, we have the code first establish arrays with 
integers 0 or 1 depending on whether the corresponding element in the 
real number array vanishes or not. These integer number arrays are 
interrogated at each integration step. Empty elements of the arrays 
are skipped; thus, a large number of useless real number multipli- 
cations are avoided. These simplifications led to a program which can 
be used in an interactive mode on a modest PC and the numbers of 
species and reactions considered can be easily changed. 
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Appendix B 


A large number of photodissociation, photoionization, ion-molecule 
reaction and recombination rates are used in this investigation. How- 
ever, only a part of these rates have a major influence on the relative 
abundances of ions and molecules discussed in Section 4. Therefore, 
we present in this paper the numerical values only for these important 
rate constants. References are given for the others. 

Photodissociation and Photoionization 

Rates corresponding to quiet solar conditions at 0.9 AU were used. 
Photoionization is of major importance for molecules that do not 
react with H 3 0 + . The rate constants for the ions listed in column 4 
of Table 2 are given in Table Bl. 

Ion-Molcule Reactions 

Rate constants were taken from the following compilations: Anicich 
and Huntress (1986), Allen et al . (1987), Prasad and Huntress (1980), 
Giguere and Huebner (1978), Huebner and Giguere (1980). The most im- 
portant rate constants, i.e. those for reactions with H 2 O and H 3 0 + , 
are given in Table 1. Atomic sulphur was introduced as a possible 
source of S + Ions (cf. text); however, reactions of ions with S were 
not included. 

Dissociative Recombination 

Dissociative recombination is usually the dominant destruction mechan- 
ism for ions that do not react with water. The corresponding rates, 
evaluated at 273 K, are given in Table B2. Rates not found in the 
literature were estimated. 
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Figure Captions 

Fl 'g- 1 « Mine mass spectra in the AMU/e range 12 to 24 obtained in 
the H-mode of the HIS sensor. Distances to the nucleus are indicated. 
The spectra obtained in the three regions of the inner coma described 
in the text are distinguished by solid lines (region I), dashed lines 
(region II), and short dashed lines (region III). 

F1 9- i - The same spectra as in Fig. 1 for the AMU/e range 25 to 35. 

Fi 9- Relative ion densities derived from the N-mode data at 1485 km 
(diamonds and bold dashed line) and 1950 km (triangles and fine dashed 
line) from the nucleus. These experimental spectra are compared with 
a model ion spectrum (crosses and solid line) which was calculated at 
the intermediate distance of 1717 km with the molecular abundances 
shown in the insert. S p indicates a point source and S £ indicates an 
extended source (R”^e"^^^) . 



Table 1 . Proton affinities (PA) and rate constants k lt k 2 



PA 

(eV) 

ki 


PA 

(eV) 

k 2 

nh 3 

8.9 

2.2 

C 2 H4 3] 

7.1 

< .001 




c 2 h 6 

6.9 

2.95 

C 2 H s 0H 

8.16 

2.8 

c 2 h 2 

6.7 

.22 

ch 3 cho 

8.07 

3.6 

HCO 

6.6 

3.0 

CH 30 H 

7.9 

2.8 

OH 

6.18 

2.89 

N 2 H 4 

7.9^ 


CO 

6.15 

2.6 

HCOOH 

7.81 

2.7 

ch 4 

5.7 

2.5 

H 2 C0 

7.5 

3.4 2 ^ 

C0 2 

5.68 

2.2 

HCN 

7.43 

4.0 

NO 

5.5 

-- 

H 2 S 

7.4 

1.9 

ch 3 

5.4 

£ .01 




n 2 

5.1 

2.5 




CN 

»v.5 

3.2 

h 2 o 

7.20 


0 2 

4.38 

< .001 




h 2 

4.38 

7.3 


Proton Affinities from Radzig and Smirnow (1980) and 
Ausloos (1974) 

Rate constants kj and k 2 , defined in the text, in 
10 9 cm 3 s 1 from Anicich and Huntress (1986) 

^from dissociation constant 
21 

J the reverse reaction has a rate of ^.23*10 9 cm 3 s 1 

3 ^kj = .06xl0' 9 cm 3 s ' 1 



Table 2. Inner coma: Ions in the 25-35 mass ranae 1 ^ 





Ions in Relation 

to Abundances of 

Parent Molecules 


Parent 

Molecules 

Major Ions 
Produced by 

31 

Long-lived J Ions 
Produced by 

Minor Ions 
Destroyed by 



Photo „i 
1 i fetime^-* 
(s) 

H 3 0 + 

Photons 

h 2 

: 0 

h 2 s 

(34) 

3.2xl0 3 

H 3 S + (35) 


h 2 s + 

(34) 






HS + 

(33) 

CH 3 OH (32) 

3.1*10 3 

CH 3 0H 2 + (33) 


CH 3 0H + 

(32) 

S 

(32) 

9 xio 5 


S + (32) 



h 2 co 

(30) 

2.9xl0 3 

CH 3 0 + (31) 


h 2 co + 

(30) 

NO 

(30) 

2.4xl0 5 


N0 + (30) 



c 2 h 6 

(30) 

4.9xl0 4 



c 2 h 6 + 

(30) 






c 2 h 7 + 

(31) 

c 2 h 4 

(28) 

1.9xl0 4 


C 2 H 4 + (28) 6] 

C 2 H 5 + 

(29) 6 ^ 

CO 

(28) 

1.3xl0 6 



co + 

(28) 






HC0 + 

(29) 

n 2 

(28) 

8 . 1 x 10 s 



n 2 + 

(28) 






n 2 h + 

(29) 

HCN 

(27) 

7.7xl0 5 

H 2 CN + (28) 


HCN + 

(27) 






CN + 

(26) 

C 2 H 2 

(26) 

2 . 6 xl 0 4 


C 2 H + (25) 4 ^ 
C 2 H 2 + ( 26 ) 5 * 6 

C 2 H 3 + 

(27) 6] 


^ Mass numbers in parenthesis 

2 -l 0.9 AU; data from Levine (1985); Kim and A'Hearn (1990) 

^ Not destroyed by H 2 O 

^ Possibly also produced by He + near and outside the contact surface 
^ Slowly reacting with H 2 0 (cf. Table 1) 

^ Various hydrocarbons may contribute to this ion 



Table 3. Molecules and Ions Included in the Coma 
Model Calculations 


Molecules 

H, H 2 , 0, OH, H 2 0 

C, C H g. , C 2 H 2 , C 2 H 4 , C 2 H 6 

N, NH 3 , HCN, N 2 , NO 

co, h 2 co, ch 3 oh, co 2 
s, h 2 s 

Ions 

H + , H 2 + , 0 + , 0H + , H 2 0 + , H 3 0 + 

C + , CH 2 + , CH 3 + , CH 4 + , CH 5 + 

c 2 h + , c 2 h 2 + , c 2 h 3 + , c 2 h 4 + , c 2 h 5 + , c 2 h 6 + , c 2 h 7 

N + , NH + , NH 2 + , NH 3 + , NH 4 + 

HCN + , H 2 CN + , N 2 + , N 2 H + , N0 + 

C0 + , HC0 + , H 2 C0 + , H 3 C0 + , CH 3 0H + , CH 3 0H 2 + 

C0 2 + , HC0 2 + 

S + , HS + , H 2 S + , H 3 S + 



Table 4. Relative contributions of ions to some important mass numbers. 
These percentages refer to the theoretical ion mass spectrum 
shown in Fig. 3. 


Mass-28 


Mass-30 

Mass-31 

Mass-32 

Mass-33 

H 2 CN + = 

69% 

N0 + = 87% 

H 3 C0 + = 100% 

CH 3 0H + = 62% 

CH 3 0H 2 + - 99% 

C0 + = 

20% 

H a C0 + = 12% 


S + = 38% 

SH + = 1% 

C 2 H 4 + = 

11% 

C 2 H 6 + = 1% 





Table Bl. Photoionization rates at 0.9 AU 


Molecule 

Ion 

Rate 

10" 7 s" 1 

H 2 0 

H 2 0 + 

4.1 

CO 

C0 + 

3.8 

NO 

N0 + 

16.0 

c 2 h 2 

C 2 H 2 + 

9.6 

c 2 h 2 

c 2 h + 

.9 

C 2 H 4 

c 2 h 4 + 

7.2 

c 2 h 4 

c 2 h 2 + 

2.5 

s 

s + 

12.0 


Data are for quiet solar condition 
(Levine, 1985; Appendix I) at 0.9 AU 



Table B2. Dissociative recombination rates 
with electrons at 273 K 


Ion 

k 

(10 _7 

cm 3 s" 1 ) 

nh 4 + 

15.9 

(a) 

h 3 o + 

10 

(a,b) 

c 2 h + 

5.7 

(a) 

c 2 h 2 + 

5.7 

(a) 

c 2 h 4 + 

10 

(d) 

h 2 cn + 

3.5 

(c) 

N0 + 

4.6 

(a) 

ch 3 0 + 

10 

(d) 

ch 3 oh 2 + 

8.8 

(c) 

h 3 s + 

3.7 

(c) 


(a) McGowan and Mitchell (1984) 

(b) Heppner et al . (1976) 

(c) Adams and Smith (1988) 

(d) estimate 
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